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Abstract 

The role of chiral transformations in effective theories modeling Quantum Chromo Dy- 
namics is reviewed. In the context of the Nambu-Jona-Lasinio model the hidden gauge 
and massive Yang-Mills approaches to vector mesons are demonstrated to be linked by a 
special chiral transformation which removes the chiral field from the scalar-pseudoscalar 
sector. The role of this transformation in the presence of a topologically non-trivial chiral 
field is illuminated. The fermion determinant for such a field configuration is evaluated 
by summing the discretized eigenvalues of the Dirac Hamiltonian. This discretization 
is accomplished by demanding certain boundary conditions on the quark fields leaving 
a finite model space. The properties of two sets of boundary conditions are compared. 
When the topologically non-trivial chiral transformation is applied to the meson fields 
the associated transformation of the boundary conditions is shown to be indispensable. 
A constructive procedure for transforming the boundary conditions is developed. 
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1 Introduction 



Since a solution to Quantum Chromo Dynamics (QCD) is not yet available one has 
to recede on models in order to explore processes described by the strong interaction. 
These models are usually constructed under the requirement that the symmetries of the 
underlying theory, i.e. QCD are maintained. In this context chiral symmetry and its 
spontaneous breaking play a key role. In this article we will explore a special chiral 
transformation when topologically non-trivial meson field configurations like solitons are 
involved. To begin with, let us briefly review the relevance of chiral symmetry on the one 
side and solitonic field configurations on the other in the context of strong interactions. 

QCD can be extended from SU(3) to SU(Nc) where Nc denotes the number of 
color degrees of freedom. It was observed by 't Hooft[[I[] that in the limit Nq — > oo 
QCD is equivalent to an effective theory of weakly interacting mesons. Subsequently 
Witten0 conjectured that baryons emerge as solitons of the meson fields within this 
effective theory. Stimulated by Witten's conjecture much interest has been devoted to 
the description of baryons as chiral solitons during the past decade @, ||. In the soliton 
description of baryons the chiral field and in particular its topological structure play 
a key role. The topological character of the chiral field especially endows the soliton 
with baryonic properties like baryon charge and spin||. This comes about via the chiral 
anomaly || which is a unique feature of all quantum field theories where fermions live in 
a gauge group and couple to a chiral field. The fermionic part of such a theory has the 
generic structure 



Z F [$\ = J D^Dty exp 



d 4 x * (ip - m ) * 



= Det (ip - rh ) (1) 

where rho denotes the current quark mass matrix which will be ignored in the ongoing 
discussions. Furthermore 

ip=i$-<$> = i($ + f)- MP R - M^P L (2) 

represents the Dirac-operator of the fermions in the external Bose-field $. $ in general 
contains vector V^, axial vector fields^ as well as scalar S and pseudo-scalar fields P 

T, = V, + A^, M = S + iP = &E£ R . (3) 

Here Pr/l = |(1 ± 7s) are the chiral projectors. Accordingly one defines left(L)- and 
right (R)-handed quark fields: ^l,r = Pl,r^ ■ The chiral field U is defined via the polar 
decomposition of the meson fields 

U = & R = exp(*6) (4) 

The chiral anomaly arises because there is no regularization scheme which simulta- 
neously preserves local vector and axial vector (chiral) symmetries. In renormalizable 
theories the chiral anomaly can be calculated in closed form and is given by the Wess- 
Zumino action 0, 

^ = Swz = ^-J [Udtff (5) 
2407T 2 Jm 5 



our notation and are anti-hermitian. 
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where the notation of alternating differential forms has been used. M 5 denotes a five 
dimensional manifold whose boundary is Minkowski space. Obviously the chiral anomaly 
is tightly related to the chiral field since the Wess-Zumino action vanishes when the chiral 
field disappears (U =1). The chiral anomaly is, however, not merely a technical artifact 
but has well established physical consequences. In the meson sector it gives rise to the 
so-called "anomalous decay processes" like e.g. it — > 27 and uo — > 3tt. In the soliton 
sector the chiral anomaly requires for Nc = 3 the soliton to be quantized as a fermion 
and endows the soliton with half integer spin and integer baryon number!||, [|. 

For many purposes it is convenient to perform a chiral rotation of the fermions[T(| |TTJ 

q = nq with n = P L £ L + P R £ R , i.e. M',../,. O../^'/../,'- (6) 
This transformation defines a chirally rotated Dirac-operator 

MpV = ^fp^! (7) 

which acquires the form 

i~p= SlHptf = i7„ (d" + W + A^ 5 ) -£. (8) 

The chiral rotation has removed the chiral field from the scalar pseudo-scalar sector of the 
rotated Dirac operator ip. As a consequence the vector and axial vector fields become 
now chirally rotated 

%-K = Ud, + V,~A^[. (9) 

Even in the absence of vector and axial vector fields in the original Dirac operator (V^ = 
A^ = 0) the chiral rotation induces vector and axial vector fields 

V ll (V ll = A ll = 0) = v ll = ^((^ + (4(1), 

A,(V, = A, = 0) = a, = i(^^-a^i). (10) 

For the soliton description of baryons the chiral field is usually assumed to be of the 
hedgehog type 

U = exp (i6(r)r • r) (11) 

The non-trivial topological structure of this configuration is then exhibited by the bound- 
ary conditions O(0) = — rnr and ©(00) = 0. The chiral field thus represents a mapping 
from the compactified coordinate space (all points at spatial infinity are identified) to 
SU(2) flavor space, i.e. S 3 — > S 3 . The associated homotopy group, n 3 (5' 3 ), is isomorphic 
to Z, the group of integer numbers. The isomorphism is given by the winding number 
(6(0) — 6(00)) /ir — —n. Assuming the unitary gauge {£} L = £r) the induced vector field 
is of the Wu-Yang form[pT[] 

v = 0, v i = iv£, vt = ^ ha f k — (12) 
2 r 

2 For this proof it is mandatory to consider flavor SU(3). 
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with the profile function G(r) given by the chiral angle 0(r) 



G(r) = -2sin 2 ^. (13) 

For odd n the topological non-trivial character of the chiral rotation is also reflected by 
a non- vanishing value of the induced vector field at r = 0. The induced axial vector field 
a, = m"r a /2 becomes 

a ~~ ( r^u \ sin6(r)\ sin6(r) 
< = nr a I e'(r) ^ J + 5 ia ^ (14) 

where the prime indicates the derivative with respect to the argument. 

The use of the chirally rotated fermions is advantageous since the rotated Dirac- 
operator does no longer contain a chiral field and its determinant is hence anomaly free. 
The chiral anomaly, however, has not been lost by the chiral rotation but is hidden in 
the integration measure over the fermion fields. In fact, as observed by Fujikawa ||12||, a 



chiral rotation of the fermion fields gives rise to a non-trivial Jacobian of the integration 
measure 

D^D^l = J(U)D^D^ (15) 
which is precisely given by the chiral anomaly 

J(U) = exp (iA) (16) 

Thus we have the relation 

-iTr\ogip= -iTr\ogip+ A (17) 

In many cases it is convenient to work with the chirally rotated fermion fields because of 
the absence of the anomaly from the fermion determinant. 

As already mentioned above the chiral anomaly or equivalently the non-trivial Jaco- 
bian in the fermionic integration measure arise due to the need for regularization, which 
introduces a finite cut-off. In the regularized theory the chiral anomaly can be evaluated 
in a gradient expansion. In leading order the anomaly is then given by the Wess-Zumino 
action (]5|). There are, however, higher order terms which are suppressed by inverse pow- 
ers of the cut-off. In renormalizable theories where the cut-off goes to infinity, these 
higher order terms disappear and the chiral anomaly is known in closed form. In non- 
renormalizable effective theories, however, the cut-off of the regularization scheme has 
to be kept finite and acquires a physical meaning, indicating the range of validity of the 
effective theory. In this case the higher order terms of the gradient expansion do not 
disappear but contribute to the anomaly which is then no longer available in closed form. 

Furthermore, when the soliton sector of such effective mesonic theories is studied it is 
not sufficient to only consider the leading and sub-leading contributions from the gradient 
expansion and one has to perform a full non-perturbative evaluation of the fermion 



determinant [13fl . The non-perturbative calculations have to be performed numerical- 



ly!^], with the continuous space being discretized. Also in the non-perturbative studies 



of the soliton sector of the effective theory the use of the chiral rotation is in many cases 



advantageous[ll]. As noticed above for the soliton description of baryons the topological 
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nature of the chiral field is crucial. Actually, topology is a property of continuous spaces 
(manifolds) and it is a priori not clear whether the chiral rotation with a topologically 
non-trivial chiral field can be represented in a finite dimensional and discretized model 
space used in the numerical calculations. In this respect let us recall that a single point 
defect in a manifold changes its topological properties already drastically. 

The present paper is devoted to a study of the chiral rotation in non-perturbative 
soliton calculations where the fermion determinant has to be numerically evaluated in 
the background of a topologically non-trivial chiral field. For definiteness we shall use 
the Nambu-Jona-Lasinio (NJL)]nj model as a microscopic fermion theory which shares 
all the relevant properties of chiral dynamics with QCD. Its bosonized version gives a 
quite satisfactory description of mesons [ ID | and also of baryons when the soliton picture 
is assumed ||16||. The organization of the paper is as follows: After these introductory 
remarks we review the importance of the local chiral rotation (|6|) for the extraction of 
meson properties from the NJL action. In section 3 we discuss the soliton solution to 
the NJL model of pseudoscalar fields with special emphasis on the choice of boundary 
conditions in the finite model space. Section 4 finally is devoted to the study of the local 
chiral rotation for topologically non-trivial field configurations and its influence on the 
soliton solution. Concluding remarks are given in section 5. 



2 Properties of local chiral transformations 

As pointed out in the introduction the chirally rotated formulation of the NJL model 
is suited to investigate properties of (axial-) vector mesons. In the present section we 
will therefore briefly review the results and illuminate the connection with the hidden 
gauge symmetry (HGS) and massive Yang-Mills (MYM) approaches for the description 
of (axial-) vector mesons. Most of the results reported in this section are taken from 
earlier works [lil 
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Nevertheless we repeat these results here in order to put our 
work into perspective and have the paper self-contained. 

In the chirally rotated formulation the bosonized version of the NJL model action 
reads 



-4njl = -ii:i\ogip+ A 



tr £ 



m~ 



4G> 



-tr 



(^-^) 2 + (4-<v)' 



with the constituent quark mass m being the vacuum expectation value of the scalar 
field £ i.e. (£) = m. Again we have discarded terms proportional to the current quark 
mass matrix. The chirally transformed (axial-) vector fields are defined in eqns @ and 
(|T0|). Next we have to face the fact that the functional trace of the logarithm in ( |T8| ) is 
ultra-violet divergent and thus needs regularization. This is achieved by first continuing 
to Euclidean space (x = — zx 4 ) and then representing the real part of the Euclidean 
action by a parameter integral 



^Trlog (f} E ~p E 



ds 

— exp 

1/A 2 S 



•$ePe 



(19) 



which introduces the cut-off A. This substitution is an identity up to an irrelevant 
constant for A — ► oo. The Euclidean Dirac operator p E is obtained by analytical con- 
tinuation of p to Euclidean space. The prescription ([19]) is known as the proper-time 
regularization IT^]. For the purpose of the present paper it is sufficient to only consider 
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the normal parts of the action. We may therefore neglect the imaginary part as well as 
the anomaly A. Thus the actual starting point of our considerations is represented by 



A 



NJL 



1/A2 T Tr ex p (- s ~APe 

1 

tr 



4G 2 



Ad 

2 



-tr E 



m 



+ (Af, - a„ 



(20) 



In order to extract information about the properties of the (axial-) vector mesons com- 
monly a (covariant) derivative expansion of the fermion determinant is performed. We 
choose to consider a covariant derivative expansion since, in contrast to an on-shell de- 



termination of the parameters [20], it preserves gauge invariance and does not lead to 
artificial mass terms. Furthermore, this procedure leaves the extraction of the axial- 
vector meson mass unique. Continuing back to Minkowski space and substituting the 



scalar field S by its vacuum expectation value yields the leading terms JD], [IT 

6m 2 



NJL 



-tr 



9v 



2 trAj 



1 

AG 2 



-tr 



V, 



A, 



'21) 



Here and A^ denote the field strength tensors 



V 



A 



dfiV u - d v % + [%, V v ] + [A„ A u ), 
d^-d^+lV^AJ + iA^K) 



(22) 



of vector and axial-vector fields, respectively. Obviously in our convention these fields 
contain the coupling constant gy which in the proper-time regularization is given by 



9v = 4tt 



2N, 



C-. 




(23) 



For the description of the pion fields we adopt the unitary gauge for the chiral fields: £ = 
£l = £r- The pions come into the game by the non-linear realization £ = exp (it ■ tz/ f). 
Then the last term in eqn (ETJ) contains the axial-yector pion mixing which is eliminated 
by a corresponding shift in the axial field: A^ — > A'^ = A^ + (ig v m 2 jYlfG^d^T ■ ir. This 
shift obviously provides an additional kinetic term for the pions and thus effects the pion 
decay constant 



f 

J TV 



with the (axial-) vector meson masses 



Mi 



My 



AM\G 2 



(24) 



M 2 v = ^r and 



4G, 



M\ = M V + 6m 2 . 



(25) 



This brief summary of known results has demonstrated the usefulness of the chiral ro- 
tation (0) especially in the context of the derivative expansion since it eliminates the 
derivative of the chiral field from the fermion determinant (^). 

The Lagrangian of the hidden gauge approach can be obtained from (^JJ) by the 
following approximation. One neglects the kinetic parts for the axial-vector field A'^ 
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which leaves this field only as an auxiliary field. This allows to employ the corresponding 



equation of motion to eliminate A'^ resulting in 



1 



2 9 f V - 



afi tr V, 



2 1 

- ~ A L tr % . 



(26) 



In the work of Bando et al.[2~I| a was left as an undetermined parameter. Here it is fixed 
in terms of physical quantities 



M\-Ml 



(27) 



Assuming the constituent quark mass m = My / Vd not only yields the Weinberg relation 
Ma = V2M V [^ but also the KSRF relation q = 2|23|. 

Alternatively one might apply the same manipulations to the formulation in terms of 
the unrotated fields @. Then the chiral field still appears in the fermion determinant 
and one has to deal with the covariant derivative 



V fl U = d,U+[V„U]-{A fl ,U}. 
The leading terms in the Lagrangian can readily be obtained [TO] 

3m -tr {vjJWrf) + ±.tr (VJ, + A 9 1 



C 



2<7 



[ 2 



AG: 



^[VZ + Al 



(28) 



(29) 



which exactly represent the massive Yang-Mills Lagrangian] 
(axial-) vector fields according to © and noting that |[T7| 



tr A 2 , = — tr (V U UV^U ] 



24]. Transforming the 



(30) 



one immediately observes that fl29"D and (|2l| ) describe the same physics. In the language of 
the NJL model the identity of the HGS and MYM approaches stems form the invariance 
of the module of the fermion determinant under the special chiral rotation (||). 

Let us next explore the behavior of the fields under flavor rotations gi, gn- These 
are defined for the unrotated left- and right-handed quark fields 



and 



(31) 



The term which describes the coupling of the quarks to the scalar and pseudoscalar 
mesons is left invariant by demanding 



- 9l^r9r 



which introduces the hidden gauge transformation h\T0 
-> h^ L g{ , 0? h^ R g R and 



(32) 



(33) 



Obviously the scalar fields transform homogeneously under the hidden gauge transfor- 
mation. In this context it is important to note that h may not be chosen independently 
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but rather depends on the gauge adopted for the chiral fields. Consider e.g. the unitary 
gauge £' L = = £. This requires the transformation property 

£ - <7i#i = h^gl (34) 

For vector type transformations gi = 9r = 9v this equation is obviously solved by h = g v . 
Contrary, for axial type transformations gi = g\i = gA h is obtained as the solution to 
QA^h = h^gA which depends on the field configuration £. Thus even for global flavor 
transformations g^ y the hidden gauge transformation h may be coordinate-dependent 
for coordinate dependent £ (x) [p5| , p6| . The unrotated (axial-) vector fields transform 



inhomogeneously under the flavor rotations 

Vp + An^gnfin + Vn + AJgl and V» - A„ -> g L (d M + - A„) g{. (35) 

It is then straightforward to verify that the flavor transformation of the rotated fields 
only involves the hidden symmetry transformation h\TQ 



*L,R^rf* L)R , %^h) (d^ + V^h and A^tfA^h. (36) 

The fact that A^ transforms homogeneously has the important consequence that one 
can build models by setting A^ = without breaking chiral symmetry. These models 
then do no longer contain axial-vector degrees of freedom. Such models have frequently 
been utilized in the investigation of p and u meson physicspB], |27|. They have also been 
successfully applied to describe baryons as solitonspHL p9], jSOfl. 



3 The NJL soliton 

In the baryon number one sector the NJL model has the celebrated feature to pos- 
sess localized static solutions with finite energy, i.e. solitons|13|, [IJJ. Here we wish to 



briefly review this solution and discuss different boundary conditions for the quark wave- 
functions. 

For static field configurations it is convenient to introduce a Dirac Hamiltonian 

H = at • p + (3 (p R mZ + Pl^ (£>£ f ) (37) 

where we have assumed the unitary gauge (i.e. Cl — ^R — 0- This Hamiltonian enters 
the Euclidean Dirac operator via 

iPfe = -d T - H. (38) 

For static mesonic background fields the fermion determinant can conveniently be ex- 
pressed in terms of the eigenvalues e M of the Dirac Hamiltonian 

= eA- (39) 

These eigenvalues obviously are functionals of the mesonic background fields. Depend- 
ing of the specific boundary condition (which fixes the quantum reference state) to the 
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fermion fields in the functional integral flU), the fermion determinant ([!]) contains in 
general besides a vacuum part A also a valence quark part ^4 val |l3 



A = A" + A val . (40) 

The valence quark part arising from the explicit occupation of the valence quark levels 
is given by 

A™ 1 = -E V ^]T , £ val [£] = N c £ (41) 

Here rj^ = 0, 1 denote the occupation numbers of the valence (anti-) quark states. The 
vacuum part is conveniently evaluated for infinite Euclidean times (T — > oo) which fixes 
the vacuum state as the quantum reference state (no valence quark orbit occupied). For 
the present considerations it will be sufficient to evaluate the real vacuum part 

A% = ^Tr \ogp E p E . (42) 

Since for static configurations one has [d T ,TC]—0 and thus ]$ e Pe = —d 2 + 7i 2 . Then it 
is straightforward to evaluate the real part of the fermion determinant in proper-time 
regularizationl 



2 Jl/A 2 s 

N C f°° d Zsr [°° ds / r 2 2] x 

The temporal part of the trace has become the z integration. As this integral is Gaussian 
it can readily be carried out yielding 

This expression allows to read off the static energy functional E[£\ since A° R ~~ > ~ TE°[£\ 
as T — > oo 

^ = ^T *£exp (-«?). (45) 



A' 



The occupation numbers 77^ of the valence quark orbits have to adjusted such the baryon 
number 

B = E " \^ M) (46) 

equals unity. The total energy functional is finally given by 

E[£] = £ val [£] + £ vac [£] - E°[£ = 1] (47) 



3 The imaginary part does not contribute for the field configurations under consideration. 
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which is normalized to the energy of the vacuum configuration £ = 1. In the chiral limit 
(m^ = 0), which we have adopted here, the meson part of the action does not contribute 
to the soliton energy. The chiral soliton is the £ configuration which minimizes E[£] and 
the minimal E[£\ is then identified as the soliton mass. 

To be specific we employ the hedgehog ansatz for the chiral field 

£(r) = exp ( % -r • f &(r)) (48) 



2 

while the scalar fields are constrained to the chiral circle, i.e. (E) = m. Substituting this 
ansatz into the Dirac Hamiltonian (|37| ) yields 



H = a • p + (5m (cosO(r) + ry 5 r • f sinO(r)) . (49) 

The stationary condition 5E[£\/5£ = is made explicit by functionally differentiating 
the energy-eigenvalues e M with respect to 



5e M 

— m 



J dtt #],(r)/3 (-sinG(r) + i 75 r ■ f cos6(r)) ^(r). (50) 



5Q(r) 

This leads to the equation of motion | 

cos6(r) tr / dVt ps(r,r)i'y 5 T ■ r = sin©(r) tr / dVt ps(r,r) (51) 

where the traces are over flavor and Dirac indices only. According to the sum ( fTTD the 
scalar quark density matrix ps(x, y) = (g(x)g(y)) is decomposed into valence quark and 
Dirac sea parts: 

ps{x,y) = p v s a \x,y)+p™(x,y) 



Technically the discretized eigenvalues e M of the Dirac Hamiltonian Ti. (|37|j39|) are 
obtained by restricting the space R3 to a spherical cavity of radius D and demanding 
certain boundary conditions at r = D. Eventually the continuum limit D — > 00 has to be 
considered. In order to discuss pertinent boundary conditions it is necessary to describe 
the structure of the eigenstates of TL. Due to the special form of the hedgehog ansatz the 
Dirac Hamiltonian commutes with the grand spin operator 

G = J + - = Z + - + - (53) 
2 2 2 y 1 

where J labels the total spin and I the orbital angular momentum, r/2 and a/2 de- 
note the isospin and spin operators, respectively. The eigenstates of 7i are then as well 
eigenstates of G. The latter are constructed by first coupling spin and orbital angular 
momentum to the total spin which is subsequently coupled with the isospin to the grand 
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spin^T]] . The resulting states are denoted by \ljGM) with M being the projection of G. 
These states obey the selection rules 




(54) 



The Dirac Hamiltonian furthermore commutes with the parity operator. Thus the eigen- 
states of 7i with different parity and/or grand spin decouple. The coordinate space 
representation of the eigenstates is finally given by 



(55) 



tfO,+) = ( igf^\r)\GG + \GM) \ ( igf^{r)\GG-\GM) \ 
~ \ /f ' +;1) (r)|G + 1G + \GM) ) + ^ -/f ^(r)\G - 1G - \GM) ) 

*«?,-) _ (ig^K^G+lG+lGM^ ( ig^\r)\G - 1G - §GM>\ 
" " ^ -f^-' l Kr)\GG+\GM) ) + \ fiG^)(r)\GG-\GM) )' [M) 

The second superscript labels the intrinsic parity Ili ntr which enters the parity eigenvalue 
via I! = (— 1) G x Ilintr. In the absence of the soliton (i.e. O = 0) the radial functions 
(r) , /^ G,+;1 ^(r), etc. are given by spherical Bessel functions. E.g. 

g f,+V( r ) = N k ^l + m/E j G (kr), /f - +;1) (r) = N k sign(E) ^1 - m/E j G+ i(kr) (57) 

and all other radial functions vanishing represents a solution to 7i(Q = 0) with the energy 
eigenvalues E = ±V^ 2 + rn 2 and parity (— 1) G . iV& is a normalization constant. 

Two distinct sets of boundary conditions have been considered in the literature. Orig- 
inally Kahana and Ripka|31|] proposed to discretize the momenta by enforcing those com- 
ponents of the Dirac spinors to vanish at the boundary which possess identical grand spin 
and orbital angular momentum, i.e. 
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This boundary condition has the advantage that for a given grand spin channel G only 
one set of basis momenta {k n c\ is involved. These k n Q make the G th Bessel function 
vanish at the boundary (jc^ncD) = 0). However, this boundary condition has (among 
others) the disadvantage that the matrix elements of flavor generators, like T3 are not 
diagonal in momentum space. If the matrix elements of the flavor generators are not 
diagonal in the momenta a finite moment of inertia will result even in the absence of a 



chiral field[I3|. Stated otherwise, in this case the boundary conditions violate the flavor 
symmetry. This problem can be cured |32| by changing the boundary conditions for the 
states with n intr = — 1 

gf^\D) = gf^\D) = gf>~*\D) = gf^\D) = (59) 

i.e. the upper components of the Dirac spinors always vanish at the boundary. The 
diagonalization of the Dirac Hamiltonian (|37D with the condition (|59|) is technically less 
feasible since it involves three sets of basis momenta {k n c-i}, {k n c+i} and {k n c} for a 
given grand spin channel. In table 1 we compare some properties of the two boundary 
conditions (|58|) and (|59D in the case when no soliton is present. The first four quantities 
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Table 1. Properties of the two boundary conditions ( pq ) and (£^) in the baryon number zero sector. 
f(r) represents an arbitrary radial function. 



Quantity Condition ( p8|) Condition fl59|) 



E M *^75«*r • f ^erfc 




A 


) sign(e M ) = 





yes 


yes 


E M ^«- (r x f) ^erfc( 


A 


) sign(e M ) 


= 


no 


yes 


E M ■ T^erfc ( 


A 


i 


si 


gn 


= 




yes 


yes 


E M • fr • r^erfc ( 


A 


) sign(e M ) = 





yes 


yes 


tr(/3/(r)) = 














yes 


no 


tr ( 7 . r /(r)) = 














yes 


yes 


(/x Ti \u) = for £^ ^ k 


V 










no 


yes 



appearing in that table show up in various equations of motion when e.g. also (axial-) 



vector mesons are included [|33]]. In case such a quantity is non-zero the vacuum gives a 
spurious contribution to the associated equation of motion. For an iterative solution to 
the equations of motion this spurious contribution has to be subtracted. It should be 
noted, however, that the relations listed in table 1 are all satisfied for both boundary 
conditions in the continuum limit D — > oo. 

So far the discussion of the boundary conditions has only effected the point r = D. 
In the context of the local chiral rotation it is equally important to consider the wave- 
functions at r = 0. As already mentioned the solutions to the Dirac equation (pE9| ) are 
given by spherical Bessel functions in the free case, = 0. Except of jo these vanish at 
the origin. In the case ^ we adopt the boundary conditions 0(0) = —nn thus no 
singularity appears in the Dirac Hamiltonian ( P?D at r = 0. Therefore the radial parts 
of the quark wave-functions may be expressed as linear combinations of the solutions to 
the free Dirac Hamiltonian. E.g. 



9f' +;1 \r) = J EV l *[Q)N k y/l + m/E kG j G (k kG r), 

k 



/f (r) = E ^[©^sign^yT - m/E kG j G+l {k kG r) (60) 
k 

where the eigenvectors V^[0] are obtained by diagonalizing the Dirac Hamiltonian in the 
free basis. It should be stressed that the use of the free basis is only applicable because 
the point singularity hidden in r • f has disappeared. If singularities show up for certain 
field configurations the basis for diagonalizing 7i has to be altered. This will be the 
central issue of the next section. 



4 The chirally rotated fermion determinant 

In the previous section we have demonstrated that the normalizable solutions to the 
free Dirac equation with spherical boundary conditions represent a pertinent basis for 
the diagonalization of the Dirac Hamiltonian with the soliton present. This property is 
based on the fact that the Hamiltonian (|I9|) is free of singularities. In this section we will 
explain how singularities appearing in a Dirac Hamiltonian influence the choice of basis 
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states. Let us for this purpose consider the Hamiltonian for the chirally rotated quark 
fields * = Q(Q)V (cf. eqns ©-(0) and ref.|TT|): 



H R = n(Q)nnUe) = a-p + /?m--(o--r)(r-r) (V(r) - -sinO(r 

2 V r 

-—(a- ■ r) sin 0(r) - -a • (r x r) sin 2 [ ^ ] (61) 
2r r V 2 / 

since in unitary gauge fi(O) = cos(0/2)+27 5 T-f sin(6/2). Obviously the ©-dependence 
in the Hamiltonian has been transferred to induced (axial-) vector meson fields. As ex- 
pected the rotated Hamiltonian, H,r, contains an explicit singularity in the 1/ra ■ (f x 
t) sin 2 f ^f^) term at r = 0. Additionally there are "coordinate singularities" in the 
expressions involving f . All these singularities appear because the "coordinate singular- 
ity" in fi(O) is not compensated by corresponding values of the chiral angle 0. Stated 
otherwise: the transformation Q(Q) with O(0) — O(oo) = — nir is topologically distinct 
from the unit transformation. Although f2(0) represents a unitary transformation it is 
then not astonishing that a numerical diagonalization 

H R V^ = e M * M (62) 

in the basis of the free Hamiltonian does not render the eigenvalues of the unrotated 
Hamiltonian, TC, (i.e. e M ^ e M ) despite the relevant matrix elements being finite. This 
finiteness is merely due to the r 2 factor in the volume element. One might suspect 
that the Hamiltonian TC2R = Vt(2Q)TiVft (20) obtained by a 20 rotation has the same 
spectrum as TC, since H.2R is free of singularities. Although this behavior is exhibited by 
the numerical solution for the low-lying energy eigenvalues, the topological character of 
the transformation has drastic consequences for the states at the lower and upper ends 
of the spectrum in momentum space. Adopting the same basis states for diagonalizing 
TC and TC2R one observes that the eigenvalues of TC2R are shifted against those of TC, i.e. 
the most negative energy eigenvalue is missing while an additional one has popped up 
at the upper end of the spectrum. Up to numerical uncertainties the eigenvalues in the 
intermediate region agree for both 7i and H,2r. This behavior is sketched in figure 1 
and repeats itself for TC^r = fi(4©)7Y^(4©). Thus the chiral rotation represents another 
example of the so-called "infinite hotel story" [34| which is an interesting feature reflecting 
the topological character of this transformation. 

Let us now return to the problem of diagonalizing Hr and restrict ourselves for the 
moment to the channel G u = + . At r = the chiral rotation 

ft(r = 0) = -z(r • f)7 5 (63) 

obviously exchanges upper and lower components of Dirac spinorsi. The corresponding 
wave-functions are given by 

n(r - 0)^+)(r - 0) - (-f^ T = °)I°! 00 M (64) 

which obviously cannot be represented by the free basis since g^ ,+,1 \r = 0) 7^ 0, in 
general. However, the lower components of the eigenstates of the free Hamiltonian in the 



J In the standard representation 75 = ( ^ q 
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TL-2R 7~L 



Fig. 1. A schematic plot of the spectrum of the rotated Hamiltonian 7Y n # = Q,(nQ) / HQ,^ (nQ) . 



G u = + channel have this property (c/. eqns (|55| , [57| ) ) . We therefore apply the rotation 
(|63|) to the eigenstates of the free Hamiltonian as well, resulting in a new basis given by 




This is, of course, no longer a solution to the free Dirac equation. Moreover, a single 
component of this spinor does not even solve the free Klein Gordan equation. At r = D 
the chiral rotation equals unity. Thus we demand the discretization condition ji(kD) = 
according to eqn (|58f) i.e. the momenta are taken from the set {k n {\. With this basis 
we have succeeded in eliminating the singularities at r = and keeping track of the 
boundary conditions at r = D. We are thus enabled to numerically diagonalize 1~Lr. 
When comparing with the eigenvalues of 7i we again encounter a form of the "infinite 
hotel story" : one state is missing in the negative part of the spectrum while an additional 
shows up in the positive part. The missing state turns out to be the one at the upper end 
of the negative Dirac sea, i.e. E « —m. Then it is important to note that in addition to 
the states with finite k the basis fl6"5| ) together with the boundary condition ji(kD) = 
also allows for the "constant state" with k = 0. In the continuum limit (D — > oo) this 
state is absent. Including, however, this state for finite D in the process of diagonalizing 
H,r finally renders the missing state. This is not surprising since application of the inverse 
chiral rotation Q^(r = 0) onto this "constant state" leads to an eigenstate of the free Dirac 
Hamiltonian with the eigenvalue —m. It should be remarked that for the free unrotated 
problem a "constant state" with eigenvalue — m is only allowed in the G n = 1~ channel. 
Although Q(Q) does commute with the grand spin operator its topological character 
mixes various grand spin channels via the boundary conditions. 





(65) 
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Accordingly the diagonalization of "Hr in the G n = channel demands the basis 
states 

N .i S iga(E)y/l-E/mj (hr)mOO)\ 

' ' 1 'l + E/m ji(A:r)|0|00> J " 



with the boundary condition ji(kD) = in order to be compatible with the Kahana- 
RipkapTI diagonalization of 7i. The additional "constant state" needed here corresponds 
to a state with eigenvalue +m of the free unrotated HamiltonianH 

We have finally been able to diagonalize the chirally rotated Hamiltonian Tin in the 
G = sectors by very tricky means. It should also be kept in mind that there are now ad- 
ditional states at the upper (from + ) and lower (from 0~) ends of the spectrum which do 
not possess "counterstates" in the G = part of the spectrum in the unrotated problem. 
For the dynamics of the problem they are of no importance because their contribution 
to physical quantities is damped by the regularization. However, their existence reflects 
the topological character of the chiral rotation. 

In the other channels i.e. G > 1 we have unfortunately not been able to construct 
a set of basis states which rendered the eigenvalues of the unrotated Hamiltonian along 
the approach described above. In the G = sector we have already seen that a "global 
rotation" — i(r ■ r)7 5 is needed for the basis states in order to accommodate the boundary 
conditions at r = 0. Furthermore a mixture of different grand spin channels appears since 
this "global rotation" deviates from unity at r = D. This problem can be avoided by 
defining a basis in the topologically non-trivial sector via 

= fi f (0)^V> (67) 

with being the solutions to the free unrotated Hamiltonian. <fr represents an auxiliary 
radial function satisfying the boundary conditions 0(0) = — ir and (f)(D) = 0. E.g. we 
may take 

(j){r) = — 7r ^1 — — \ exp (—tmr) (68) 

with t being a free parameter. The diagonalization of 7{r in the basis is equivalent 
to diagonalizing 

a ■ p + m(3 (cos 4>{r) + rys(T ■ f ) sin 4>(r)) 

1 



0'(r) + - sin(6(r) - </>(r)) 
r 



(T • ?)(t • f) 



2 

+ — sin(0(r) - 6(r))(cr • t) [1 - cos(6(r) - 0(r))l a ■ (f X r) (69) 

2r 2r 

in the standard basis {^o}- At this point it should be stressed again that <fi is not 
a dynamical field but merely an auxiliary field which transforms the basis such as to 
eliminate the singularities from the Dirac Hamiltonian. This property is completely 
determined by the boundary values of <p. Thus the results ought to be independent of 
the parameter t. 

For the numerical investigation we have first considered the energy eigenvalues of Tin 
employing the rotated basis flB7p for the diagonalization and compared to the eigenvalues 



5 The additional "constant states" thus do not alter the trace of the Hamiltonian. 
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of 7i. In a wide range t ~ 0.1 ... 1.4 the difference for the eigenvalues which are contained 
in the interval (—5m, 5m) is negligible, i.e. less than tenth of a percent. For the higher 
grand spin channels the window for t which yields reasonable agreement lies somewhat 
higher i« 0.8. ..1.6. Summing up the energy eigenvalues according to eqn (^) in order 
to obtain the vacuum part of the soliton energy the deviation compared to the unrotated 
formulation is about 0.5MeV. This is, of course, negligibly small since the inherited mass 
scale is given by the constituent quark mass, m which is of the order of several hundred 
MeV. By using the locally transformed basis (|67|) the wave-functions corresponding to 
the eigenstates of Tin agree reasonably well with the rotated wave-functions fi(0)\l/ M °f 
the original Dirac Hamiltonian 7i. It should be noted that a large number of momentum 
states is required to numerically gain this result. This is not surprising since in order to 
represent the functional unity an infinite number of momentum states is needed. 

Before turning to the detailed discussion of the equation of motion in the rotated 
system we would like to mention that the self-consistent profile being obtained from the 
unrotated problem also minimizes the soliton mass in the chirally rotated frame. Stated 
otherwise: each change in this profile function leads to an increase of the energy obtained 
from the eigenvalues of TCr. 

As in the formulation with in the unrotated frame the equation of motion is gained 
by extremizing the energy functional fl4~7[). The energy eigenvalues in the rotated frame, 
however, exhibit a different functional dependence on the chiral field. The functional 
derivative of these eigenvalues with respect to the chiral angle reads 

/ dn + 7 cos0 ( r ))^ 2 ^( r )( <T • f )( r ■ 

-~cos9(r)f \{r)(j ■ r^ M (r) - ^sin0(r)^(r)a • (f x t)^(i-)} (70) 

with ^(r) being the eigenstates of the rotated Dirac Hamiltonian (^l|). The derivatives 
(frO|) enter the stationary condition for the energy functional resulting in the equation of 
motion 

A L (r) + A T {r)cosQ{r) - 1/(r)sin0(r) = 0. (71) 

In order to display the radial functions Al,Aj< and V it is convenient to introduce the 
charge density pc = (g(x)g(y)"!") = p™ 1 + p™ c involving the eigenstates of Hr (cf. eqn 

mm- 

pc\^y) = E^( x )^l(y) si s n ( e /J 

) ^i(y)sign(e M ). (72) 

The radial functions A^,At and V are of axial-^^r) and vector(l/) character 

A L (r) = ^tr jdnr 2 p c (r, r)(o" • f)(r • f) (73) 
A T (r) = tr Jdnpc(T,T)[(tT-r)(T-T)-(*-T)) (74) 
V(r) = tr [ dtt p c (r,r)ct ■ (f x r). (75) 



5<d(r) 



Pc C (x,y) = — E^Werfc 



A 
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It should be noted that Al does not depend on the auxiliary field <p. It is then straight- 
forward to verify that for any meson configuration Al and At satisfy the relation 



A L (r = 0) = (-l) k A T (r = 0) (76) 

where k is defined by the value of the auxiliary field <fi at the origin 0(r = 0) = kn. The 
vector type radial function V vanishes at the origin. Thus the equation of motion (71) 



together with the relation ( [76] ) yield the boundary condition 0(r = 0) = (2n + 1)tv for 
k = 1. This is stronger than the boundary condition derived from the original equation 
of motion ([51]) which also allows for even multiples of 7r for 0(r = 0) since tr / dQ ps(r = 
0, r = 0)i75T • f = 0. Assuming the Kahana-Ripka pl| boundary conditions for the 
unrotated basis states similar considerations for r = D show that G(r = D) = 21tc 
since A^ir = D) = — Ax(r = D) as long as 4>(D) = 0. Obviously the topological charge 
associated with the chiral field in the hedgehog ansatz (0(r = 0) — G(r = D))/ir can 
assume odd values only when k is odd. Especially <fi = is prohibited in the case of unit 
baryon number. Thus the study of the boundary conditions in the chirally transformed 
system corroborates the conclusion drawn from investigating the eigenvalues and -states 
of Tin that it is mandatory to also transform the basis spinors and in particular the 
boundary conditions to the topological non-trivial sector. 

Before turning to the discussion of the numerical treatment of eqn ([71]) we would like 
to make the remark that substituting the transformation ^(r) = S7^(@) x I / At (r) into the 
original equation of motion ([51]) does not result in the relation (|70|) but rather yields the 



constraint 



A 



^(r) T5 T-f^( r ). (77) 



Thus the equation of motion ([n]) cannot be obtained by transforming the states ^ M 
but only by employing the Dirac equation in the rotated frame (|62j) to extremize the 
energy functional. The constraint (77) does not represent an over-determination of the 



the system since infinitely many states are involved. 

Due the transcendent character of the equation of motion in the rotated frame (|7l| ) 
solutions cannot be obtained for arbitrary values of the radial functions Al,At and V. 
E.g. for large distances — > requires \Al\ > \At\ in order to find a solution to eqn (|7T[). 
Thus the treatment of the NJL soliton in the chirally rotated frame is not very well suited 
for an iterative procedure to find the self-consistent solution. The reason is that a small 
deviation of the radial functions Al, At and V from those corresponding to this solution 
can render eqn ( |TT| ) indissoluble for 0(r). Then it is not unexpected that at large distances 
the solution to the rotated equation of motion (|7T| ) becomes unstable and the original pro- 
file function cannot be reproduced for r > 2fm. For smaller values of r the original profile 
function is well reproduced. In figure 2 this behavior is displayed. The self-consistent 
solution to eqn ( |5T| ) serves as ingredient to evaluate the radial functions A L ,A T and V. 
The solution to eqn Q7ID is then constructed and compared to the original profile function. 
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Fig. 2. Comparison of the self-consistent profile in the unrotated formulation (dashed line) and the 
solution to eqn ( |7l] ) (solid line). 



5 Conclusions 

We have investigated the role of chiral transformations for the evaluation of fermion 
determinants. When these transformations are topologically trivial they provide a use- 
ful tool to evaluate the chiral anomaly. Furthermore they can be used to demonstrate 
the equivalence between the hidden gauge and massive Yang-Mills approaches to vec- 
tor mesons. A generalization to the case when chiral fields have a topological charge 
different from zero is not straightforward. Even though the special transformation we 
have been considering is unitary its topological character prevents the eigenvalues and 
-vectors of the original Dirac Hamiltonian to be regained from the rotated Hamiltonian 
unless the boundary conditions are transformed to the topologically non-trivial sector 
accordingly. Furthermore we have observed that the stationary conditions to the static 
energy functional in the topologically distinct sectors are not related by the transforma- 
tion of the equation of motion. The boundary conditions for the chiral field obtained 
from the stationary condition have been found being invariant under the chiral rotation 
only when the basis quark fields are taken from the topological sector associated with 
the chiral transformation. Diagonalization of the rotated Dirac Hamiltonian in this basis 
can be reformulated into a problem where the induced vector fields ([TDD belong to the 
topologically trivial sector (cf. eqn In order to diagonalize the resulting operator 

(|B5p the standard basis ^] may be employed. 

These explorations on chiral transformation properties of the fermion determinant 
when the Dirac Hamiltonian contains topologically non-trivial (axial-) vector mesons 
may prove to be very helpful when considering the chiral invariant elimination of the 
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axial-vector fields as described at the end of section 2[2f|. In that case the chirally 
rotated Dirac Hamiltonian becomes as simple as 



G(r) 

H K = a-p + /?m + — a-ffxr). (78) 

r 

Here G(r) refers to a dynamical field and should not be confused with the induced vector 
field discussed in eqn (|T3"|). The chiral field then only appears in the purely mesonic part 
of the energy functional 

E m = ^- Jdr |(G(r) + 1 - cos6(r)) 2 + K 2 Q'{r) + sin 2 6(r)| . (79) 



0(r = 0) = — 7r then implies that G(r = 0) = — 2|]2~8l |29fl . One thus has to deal with a 



Dirac Hamiltonian which contains a topologically non-trivial vector meson field. Thus 
(f78|) cannot be treated using the standard basis |3l], [32| but rather by employing techniques 
which are analogous to those developed in section 4. Investigations in this direction are 
in progress. 

Finally we would also like to remark that this kind of singularities does not only 
appear when the Dirac Hamiltonian is considered. Such topological defects have also 



caused problems when fluctuations off vector meson solitons were invest igated||30||. In 
that case the boundary conditions for the vector meson fluctuations had to undergo a 
special gauge transformation which corresponds to the transformation of the basis quark 
spinors described here (|67|). 
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